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We study biexcitonic states in two tunnel-coupled semiconductor quantum dots and show that 
such systems provide the possibility to produce polarization-entangled photons or spin-entangled 
electrons that are spatially separated at production. We distinguish between the various spin config- 
urations and calculate the low-energy biexciton spectrum using the Heitler-London approximation 
as a function of magnetic and electric fields. The oscillator strengths for the biexciton recombi- 
nation involving the sequential emission of two photons are calculated. The entanglement of the 
polarizations resulting from the spin configuration in the biexciton states is quantified as a function 
of the photon emission angles. 
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Entanglement has been identified as an essential re- 
source for many applications in the recently developed 
field of quantum communication and quantum compu- 
tation [y. Several quantum communication schemes 
have already been successfully implemented with pairs of 
polarization-entangled photons produced by parametric 
down-conversion jl], e.g. the faithful transmission of a 
quantum state (quantum teleportation), entanglement- 
assisted classical communication (e.g., quantum super- 
dense coding), and the production of a secure crypto- 
graphic key (quantum key distribution). Recently, there 
has also been growing interest in solid-state implemen- 
tations of quantum computation using the electron spin 
as the qubit j2j , as well as quantum communication with 
spin-entangled electrons. Superconductor-normal junc- 
tions in combination with quantum dots (QDs) have 
been suggested as a device for producing entangled elec- 
trons Q . Still, the efficient and deterministic production 
of both entangled photons and electrons poses a theoret- 
ical and experimental challenge. In the case of photons, 
the use of electron-hole recombination in a single QD 
was recently suggested |Q, || . Non-resonant excitation of 
a QD is expected to produce pairs of entangled photons 
with an efficiency (production rate/pump rate) that is 
about four orders of magnitude bigger than for paramet- 
ric down-conversion 1^. 

In this paper, we study the production of polarization- 
entangled photons, or, alternatively, spin-entangled elec- 
trons, using the biexcitonic ground state in two tunnel- 
coupled QDs. For this purpose we study the low-energy 
biexcitonic states in coupled QDs, determining their en- 
ergy spectrum and their optical properties. We concen- 
trate on the spin configuration of the calculated states, 
being related to the orbital wavefunction via the Fermi 
statistics which is implemented in a Heitler-London (HL) 
ansatz for electrons and for holes. As a special quality of 
a double dot, we find that in the (spin-entangled) biex- 
citonic ground state, the biexciton favors a configuration 
with each QD occupied by one exciton, thus providing a 
basis for the separation of the entangled particles. Even 
though coupled QDs are usually separated by a distance 
less than the wavelength of the emitted light, it might still 
be possible to directly detect the photons at separate lo- 



cations. It can e.g. be expected that due to anisotropics 
the two dots have different preferred emission directions 
inclosing a fixed angle. Two subsequent photons, which 
are emitted with a time delay given by the exciton life- 
time, could then be detected separately in the far field. 

In contrast to our calculations, earlier studies for quan- 
tum computation or entanglement production with exci- 
tons in QDs concentrate on single QDs |§i R |] 

and/or on charge degrees of freedom (neglecting spin) |q, 
g, 1 1, |10|, . Also, instead of a pure electrostatic inter- 
dot coupling 1^, 0, we take into account the tunneling 
of electrons and holes between the coupled QDs. 

Biexcitons consist of two bound excitons which them- 
selves are formed by a conduction-band electron and 
a valence-band hole in a semiconductor, bound to- 
gether by the attractive Coulomb interaction. Follow- 
ing the theory of excitonic absorption in single QDs , 
the biexcitonic states in single QDs have been investi- 
gated H, H, ill' S and single ex- 
citons in coupled QDs have been observed in experi- 
ment Recently, spin spectroscopy of excitons 
in QDs was performed using polarization-resolved mag- 
nctophotoluminenscence . Two regimes can be distin- 
guished in the discussion of excitons confined in QDs . 
In the weak confinement limit ax ^ ae,ah, where ax is 
the radius of the free exciton and Oe , a^ the electron and 
hole effective Bohr radii in the QD, an exciton can (as in 
the bulk material) be considered as a boson in an external 
confinement potential. In the case of strong confinement 
ax ^ ae,ah, electrons and holes are separately confined 
in the QD and the bosonic nature of the electron-hole 
pair breaks down. Since, e.g., in bulk GaAs ax ~ 10 nm, 
we are in an intermediate regime ax ~ Oej ah for typi- 
cal QD radii. Here, we start from a strong confinement 
ansatz, i.e. from independent electrons and holes (two of 
each species), and then use the HL approximation to in- 
clude the Coulomb interaction and the tunneling. Unlike 
for bulk biexcitons, where the HL approximation fails for 
some values of ^ = TOe/m/i we are here in a differ- 
ent situation — much more similar to the H2 molecule — 
because the single particle orbitals are defined by the 
strong QD confinement, the latter playing the role of the 
( "infinitely" heavy) protons of the H2 molecule. 
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We obtain the low-energy (spin-resolved) biexciton 
spectrum in which the electrons and holes each form ei- 
ther a spin singlet or triplet. Subsequently, we calculate 
the oscillator strength, being a measure for the optical 
transition rates. The spin of the biexciton states relates 
to two different states of the polarization-entangled pho- 
ton pair produced in the recombination. We quantify 
the entanglement of the photon pair depending on the 
emission directions. The variation of the spectrum and 
the oscillator strengths due to magnetic or electric fields 
allows to use ground-state biexcitons in tunnel- coupled 
QDs as a pulsed source of entangled photon pairs. 

We model the biexciton (two electrons and two holes) 
in two coupled QDs by the Hamiltonian 



(1) 



a—e,h i—1 



where hai = (Pai + qaMT^ai)/c)'^ f^rua + Vaivai) is the 
single-particle Hamiltonian for the i-th electron (a = e, 
Qe = — e) or hole (a = h, qh = +e) in two dimen- 
sions (2D) with coordinate Vai and spin S"^ The po- 
tential Va{x,y) = niaUjlKx'^ - a^fl^a^ 2/^]/2 de- 
scribes two QDs centered at {x = ia,y = 0), sepa- 
rated by a barrier of height maW^a^/S. Electrons and 
holes have effective masses nia and confinement ener- 
gies hwa ■ The Coulomb interaction is included by He = 
i^/'^)J2ia,i)=£{i3.j) qaqp/K-lYai-rfSjl, with a diclcctric con- 
stant K (for bulk GaAs, k = 13.18). A magnetic field 
B in z direction leads to orbital effects via the vector 
potential (in the symmetric gauge) A = B{—y, x,0)/2 
and to the Zeeman term Hz = J2ai9at^BBS"'^, where 
Qa is the effective g-factor of the electron (hole) and 
fiB is the Bohr magneton. Restricting ourselves to 
the low-energy physics of QDs filled with few parti- 
cles, we can assume approximately 2D parabolic con- 
finement. We assume the simultaneous confinement of 
electrons and holes which can be realized e.g. in QDs 
formed by thickness fluctuations in a quantum well ||] 
or by self-assembled QDs [|4[ A particle in a single 
QD is thus described by the Fock-Darwin (FD) Hamil- 
tonian h^"'{rai) ||26|| , comprising a harmonic potential 
iaUj'^[{x^a)'^ +y'^]/2 and a perpendicular mag- 



'(r) 

netic field. In prospect of the HL ansatz below we write 
the single-particle part of the Hamiltonian Eq. (|l|) as 
'EJK''{^al)+h+''{r^2)]+Hw{{rc.^}) = Ho+Hw, where 



Hw{{ra^}) = EaE,^a(ra.)-«r(ral) 



(ra2)]. 



An in-plane electric field E = ey is described by He = 
e £ (j/ei + J/e2 — J/ft,i — 2/ft,2) and can be included in Hf). We 
put e = here and discuss the case e ^ below. 

The valence band is assumed to be split into well- 
separated heavy and light hole bands and only heavy- 
hole excitations are considered in the following. The FD 
ground states \D)a in the QD D = 1,2 which are used 
to make a variational HL ansatz are |26| 



{{x±af 



2. iqaay 



(2) 



where the upper (lower) sign holds for D = 1 (2), Ib = 
y^hc/eB and ba — a/I -|- {eB /2cmaUJaY . 

We now make a strong confinement ansatz by con- 
structing two-particle orbital wave functions for electrons 
and for holes separately according to the HL method, i.e. 
a symmetric (|s)" = |/ — 0)", spin singlet) and an anti- 
symmetric (1^)" = \I = 1)", spin triplet) linear combina- 
tion of two-particle states \DD')a = \D)a ® 



|/)" = A^a/(|12)« + (-l)^|21)„), 



(3) 



where N^i = 1/^2(1 + (-1)^|5„|2) and So. = a(l|2)a 
denotes the overlap (or tunneling amplitude) between the 
two orbital wave functions |1)q, and \2)a. We continue by 
forming the four biexciton states \IJ) = \I)'^^\J)'^, where 
/ = (1) for the electron singlet (triplet) and J = (1) 
for the hole singlet (triplet). The energies 



Eij = (IJ\H\IJ) =E" 



E?j, 



(4) 



with Efj = {IJ\Ha\IJ), can be calculated analytically. 



In units of tvjjp , we find Eq = E^ 



Efj 
1 



J J - 2{be + bh/rj), where 
{^lBB/hwe)Y.ai9aST, and 



1 



4 



1 



(5) 



where 2d = 2a/ a^. is the dimensionless inter-dot distance, 
tte = ^/KfmJjJe is the electronic Bohr radius, Se ~ 
exp {-<P [2b, - I /be]), Sh = exp {-dP [2bh - l/bu] /^v), 



Nij = N'^jN'^j, and ^=me/mh. For Efj, we find 



Efj = 



Eee + i-lfEee , E, 



hh 



{-lYEhh 



Ex + Eeh + {-^'^Y SeExe 

+ {-lYShExh + {~lY+' SeShExeh 
where we have used the abbreviations 



i + i-iysi 

+8Nij ' 



(6) 



Ex 

Eeh 
Ex a 
Exeh 



Xae-lip{-baCp/Xa) Iq (bad'^/Xa) , (7) 



'bry f bad'^ 
-5a exp )/o 



1 

ba 



-cVb, 

Ex exp (-^2) /o (M^) , 
2SaExexp{~bdy4bl) h {bd^l^bl) 
SeShEx {exp (6idV2)/o (6idV2) 
+ exp (MV2)/o(MV2)}- 



,(8) 

(9) 
(10) 

(11) 
(12) 



Here, Iq{x) is the zeroth-order modified Bessel function, 
c = e^-\/7r/2 /Kajujjf, is a dimensionless parameter char- 
acterizing the Coulomb interaction, Xe = Xu = Cv, 
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b = 2bebh/{bhHr]be), bi = be-l/be+[bh - 1/bh] /^r], and 
h - {[be - l/be] [bh - 2^776,] +6, [bn - l/b^d/ib^+Cvbe)- 
Fig. |l| shows the biexciton energies Ejj {I , J — 0,1 — s,t) 
in the double QD as a function of an apphed external 
magnetic field in z direction. The Zeeman interaction 
Hz causes an additional level splitting of w 0.02 hoje/T 
(assuming \ge\ ~ \gh\ ~ 1) for the triplet states with 
J2i 7^ which is not shown in Fig. ^ The electron- 
hole exchange interaction for the GaAs QDs considered 
here is reported to be only on the order of tens of ^eV | |27[ | 
and can therefore be neglected. The self-consistency of 
omitting excited single-QD states in the HL ansatz can 
be checked by comparing the energy Efj + EYj to the 
single-QD level spacing. This criterion is fulfilled for 
inter-dot distances 2a > 20 nm. In addition to the 
HL states \IJ), we consider the double occupation states 
\DDDD) for which all four particles are located on the 
same QD D = 1,2. Their energies are given by i? = 
EO + E^ + + E^, with E^ = 3(1/62 ^ ^i^2^y^Q^2^ 

and E^ = ci^/K + - 4%/6). 

We proceed to the calculation of the oscillator 
strengths of biexciton-exciton and exciton-vacuum tran- 
sitions. The oscillator strength / is a measure for the 
coupling of exciton states to the electromagnetic field and 
is proportional to the optical transition rates. For a tran- 
sition between the N + 1 and N exciton states |iV + 1) 
and |iV), the oscillator strength is defined as 



fN+l,N — 2\pNii\f /mofiUJN+l,N, 



(13) 



where mo is the bare electron mass, hLON+i,N = ^'jv+i — 
En, and pjvkA = (iV -h l|ekA ■ p|A^), where ekA is the 
unit polarization vector for a photon with momentum k 
and helicity A = ±1, and p is the electron momentum 
operator. For p^rkA we find in the dipole approximation 



E/R(0, 
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FIG. 1: Biexciton energies in units of hwe for (a) r) — Ue/coh ~ 
1/2, (b) r? = 1/C = 1-67 (fle = ah), in a 2D GaAs system 
(me = 0.067mo, mhh = 0.112mo), hu>e = 3 meV, and d — 0.7. 
The plotted HL energies Eij are -Ess (solid line), -Est (short- 
dashed), -Ets (dot-dashed), and Ett (dotted), neglecting the 
Zeeman energy. The exchange splittings Etj — Esj, J — s,t, 
for electrons are larger than for holes {En — Ejs, I — s,t) in 
(a) where r;^ < 1, but of the same order in (b) (jy^ = !)• At 
B — 0, \ss) has the lowest energy, while for larger B, there is 
a crossover to a \tt} ground state. Double occupation of a QD 
(long-dashed line) becomes more favorable with increasing rj; 
in (a), E > Eij, I,J = s,t, while in (b), E is smaller than 
some of the Ejj for small B. 



Ga ^ 271 /k {aa ~ 20 nm, 27r/fc « 1 /im), 

PJVkA = [{N+l)f ^^-^(^) Jd^rl[d^nd^s, (14) 

{o-i,rj},cr i,j 

x$Ar({r,, crj; {sj, rj})$^_^i({ri, crj, r, cr; {s^, t,}, r, a}), 

where is the ./V-exciton wavefunction, depending on 
the conduction-band electron (valence-band hole) coor- 
dinates Ti (sj) and their spins (r^) = 1...N). 
The coordinate and spin of the electron and the hole 
created or annihilated during the optical transition are 
denoted by r and cr. The inter-band momentum ma- 
trix element for a cubic crystal symmetry is given by 
Max{9) = ekA-p™(a-) ^ pcv{cos{9)-aX)/2 =pcvmax{0), 
where 9 is the angle between k and the normal to the 
plane of the 2D electron system (assuming that the latter 
coincides with one of the main axes of the cubic crystal), 
and Ep = 2plJmo {= 25.7 eV for GaAs). 

According to Eq. (p^, the orbital momentum ma- 
trix element for transitions from the exciton vacuum 
|0) to an exciton state \X) = \D)e \D)h = \DD) 
in one QD (or for the optical recombination of \X)) 
is po = M,x{9) /c;3r$t(r,r) = M,x{9)Ceh- The exci- 
ton wave function is denoted by $i(re,r;i) = (re,r^|Ar). 
From this, we find for the oscillator strength 



fxfl — 



2bo| 



E, 



= ^-^M,x{9f\Cekf 



(15) 



and Ceh = 2^/£,ribebh/{bh + ^rjbe)- In Fig. ga we plot 
IC'eft P = ///o as a function of the magnetic field, where 
/o = Epm„x{9Y / Eg denotes the oscillator strength for 
(bulk) inter-band transitions, equating huJx,o with the 
band-gap energy Eg. Since we have made a strong con- 
finement ansatz, the obtained oscillator strength is inde- 
pendent of the QD volume V . For weak confinement, one 
would expect f oc V. Fig. ^ shows the suppression of 
the exciton transition rate by an electric field. 

The momentum matrix clement pi for transitions from 
an exciton state \X) to a biexciton state \XX) is given by 
Pi = -2M„x{9) Jd^red^rhd^r^*2ire,r;rh,r)'i>iire,rh). If 
the recombining electron and hole are on the same 
QD, the integral over r yields Ceh, otherwise S^h = 
Ceh exp (-2d2 [b^ - ^T^libh + ^vbe)]) ■ 



f/fn 




2 3 4 
B/T 



50 100 150 200 250 

E/(mV/um) 



FIG. 2: Oscillator strength fx,o for GaAs QDs in units of 
/o as a function of (a) the magnetic field B (in Tesla) at 
E = and (b) the electric field E (in mV//j,m) at -B = 0, 
with r) = uOe/L^h ~ 1/2 (solid line), rj — 1/^ (dashed), rj — 4 
(dotted). For f) = the B field has no effect on Ceh- 
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We give here our result for pi for a transition between 
the HL biexciton states \XX) = \IJ) with one exciton 
on each QD and a single exciton in the final state \X) — 
\DD), a single exciton on dot D = 1,2, 



\{IJ\ei,x- P\DD)\ = 2M,xiO)y^ 

Approximating hcoxx.x ~ Eg, we plot the corresponding 
oscillator strength versus B and E in Figs. |a and |b. 

Results for fxx,x, also including the (biexciton) dou- 
ble occupation state \DDDD) and various final (exciton) 
states, will be given elsewhere [|8l. 
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FIG. 3: Oscillator strengths fxx,x for transitions between 
the biexciton states = \IJ) and a single remaining ex- 

citon on one QD in units of /o as a function of (a) the magnetic 
field B (in Tesla) at i? = and (b) the electric field E (in 
mV/^m) at _B = 0. The parameters were chosen for GaAs 
with rj = uJe/iOh — 1/2. The line styles correspond to those 
for Eij in Fig. |l|. 

The main effect of an electric field is to spatially sep- 
arate the electrons from the holes which leads to 
a reduction of the oscillator strengths |^ (cf. Figs. ||b 
and ^d). Hence, the optical transition rate can be effi- 
ciently switched off and on, thus allowing the determin- 
istic emission of one photon pair. 

Transformation of a HL biexciton state |/J) into the 
basis of two coupled excitons yields a superposition of 
dark {Sz — ±2) and bright {Sz = ±1) exciton states. 



The emitted photon states are (up to normalization) 

\Xlj) (X 1 + 1,01)1-1,^2) + (-l)'+'Vl,^l)l + l,^2), 



(17) 



where |ct,6') = iV(6')(m<^^+i(6')|CT+) +to^ _i(6')|ct_)) is the 
state of a photon emitted from the recombination of an 
electron with spin Sz = cr/2 = ±1/2 and a heavy hole 
with spin Sz ~ 3cr/2 in a direction which encloses the 
angle 6 with the normal to the plane of the 2D electron 
and hole motion. The states of right and left circular 
polarization are denoted |(T±). 

The entanglement of the state (|l^) can be quantified by 
the von Neumann entropy E. For \ss) or \tt) and emission 
of the two photons enclosing an azimuthal angle (j) — or 
TT, we obtain E — log2(l -I- X1X2) — X1X2 log2(a;ia;2)/(l -I- 
X1X2), where Xi — cos^ {6i). Note that only the emission 
of both photons perpendicular to the plane {9i —92—0) 
results in maximal entanglement {E = 1) since only in 
this case \ + l,9i) is orthogonal to \ — l,9i). In particular, 
the two photons are not entangled {E — 0) if at least 
one of them is emitted in-plane {9i = tt/2). To observe 
the proposed effect, the relaxation rate to the biexciton 
ground state must exceed the biexciton recombination 
rate. That such a regime can be reached is suggested by 
experiments with low excitation densities, see e.g. [ pO[ . 
Then, an upper limit for the pair production rate is given 
by (tx +Txx)~^ , where tx,xx is the (bi)cxciton lifetime. 

Conversely, spin-entangled electrons can be produced 
by optical absorption followed by relaxation of the biex- 
citon to its ground state. After each QD has been filled 
with an exciton, the recombination can be suppressed 
by an electric field. Having removed the holes, the elec- 
tron singlet and triplet could then in principle be distin- 



guished by a subsequent interference experiment |31 
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